Non-Markovianity as a function of system size 
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We study the non-Markovianity of the dynamics of open quantum systems as a function of system size. We 
investigate two distinct quantum processes, namely dephasing and relaxation, by considering independent and 
common environments. For the zero temperature case, we present an efficient method to evaluate a measure 
of non-Markovianity recently proposed by Luo, Fu, and Song [S. Luo, S. Fu, and H. Song, Phys. Rev. A 
86, 044101 (2012)]. We perform a comprehensive analysis of non-Markovianity for the case of a single qubit, 
including the determination of the optimal initial states, that is required for the calculation of this new measure, 
as a function of the parameters of the environment. We demonstrate that, for independent environments, the 
degree of non-Markovianity for dephasing and relaxation processes can increase with the number of the qubits in 
the system. Furthermore, we discuss the same processes for two qubits interacting with a common environment. 
In this case, we show that while the amount of non-Markovianity increases as compared to the single qubit case 
when a common relaxation process is considered, our analysis suggests that the same is not generally true for 
the dephasing process. 
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I. INTRODUCTION 

The concept of non-Markovianity is a prominent aspect of 
the dynamics of open quantum systems that has been attract- 
ing both theoretical and experimental attention [1]. It has been 
shown that non-Markovianity can be used as a tool in quantum 
protocols [2], can be employed to take advantage in quantum 
metrology [3], and can be exploited in quantum key distribu- 
tion [4]. Other concepts behind the non-Markovian dynamics 
have also been investigated, for example, the influence of en- 
vironment size [5] and the possibility to pursue new quan- 
tum technologies by using non-Markovian effects [2]. Al- 
though all those efforts have been performed, measuring non- 
Markovianity is complicated and generally very small systems 
have yet been considered [6]. However, the real usefulness of 
a quantum system for computation or simulation is only ap- 
preciable in the limit of large-scale information processing. 
Therefore, it is critical to understand the scaling properties of 
non-Markovianity for multipartite systems. To the best of our 
knowledge, this is the first study concerning the scaling prob- 
lem for a non-Markovianity measure. 

Recently, various measures for quantifying the degree of 
non-Markovianity of the dynamics of an open quantum sys- 
tem [7-10] have been introduced in the literature; however, 
there is no consensus on what precisely determines the non- 
Markovianity of a dynamical quantum process. It has been 
demonstrated that the conclusions drawn from different mea- 
sures might not agree depending on the considered physi- 
cal model. One of the most widely used measures of non- 
Markovianity has been introduced by Breuer, Laine and Piilo 
(BLP) [7]. In their seminal paper, they claimed that infor- 
mation flows only from the system into the environment for 
a Markovian process and the information flow can be mea- 
sured by the trace distance of two arbitrary quantum states, 
which probes the distinguishably between them. To imple- 
ment this measure, one needs to perform an optimization by 



checking the dynamics of the trace distance for a huge num- 
ber of initial sates. Thus, this procedure is very demanding 
and almost impracticable when dealing with multipartite sys- 
tems. Moreover, Rivas, Huelga and Plenio (RHP) have con- 
structed a measure of non-Markovianity that quantifies the de- 
viation from divisibility for a dynamical map [8], which is 
also hard to be implemented in general. In fact, the evalua- 
tion of almost all known non-Markovianity measures is very 
complicated for multipartite systems. In this work, in order 
to overcome this difficulty, we present an efficient method to 
evaluate a measure of non-Markovianity, recently proposed by 
Luo, Fu, and Song (LFS), based on the non-monotonical be- 
havior of the quantum mutual information for non-Markovian 
processes [9]. This measure coincides with other important 
ones such as BLP and RHP measures for quite general cases 
and can be straightforwardly extended for studying multipar- 
tite systems. 

When a quantum system interacts with its environment, 
there are two important physical processes that must be con- 
sidered: relaxation and decoherence (here called dephasing). 
While relaxation is associated to a process involving loss of 
energy, dephasing is associated to the loss of purity without 
any energy exchange between the system and its surround- 
ings. Indeed, the density matrix population is constant in time 
for the latter process. In this work, we explore both processes 
considering independent and common phase damping (PD) 
and amplitude damping (AD) channels to describe dephasing 
and relaxation, respectively. For a zero temperature environ- 
ment, we link the LFS measure to the rate of change of the sys- 
tem entropy S{p''{t)) and the environment entropy S{p'^{t)). 
In fact, we show that a quantum process is non-Markovian if 
the time derivate of S{p'' (t)) is greater than the time derivative 
of S{p^ (t)) at any instant of time. We present a detailed analy- 
sis of the evaluation of the LFS measure for a single qubit, and 
discuss the behavior of the optimal initial states as a function 
of the parameters of the environments. Moreover, we demon- 
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strate that the degree of non-Markovianity can increase as a 
function of the number of qubits for the cases of independent 
PD and AD channels. We also investigate the non-Markovian 
behavior of global PD and AD channels for a system of two 
qubits. Lastly, we compare the results obtained for the LFS 
measure and its simplified version, whose evaluation does not 
include a complex optimization procedure. 

II. MEASURING NON-MARKOVIANITY 

The definition of the LFS measure of non-Markovianity [9] 
is based on the following: suppose that we have a quantum 
system in a Hilbert space H, and a quantum process A(t) gov- 
erning the dynamical evolution of the considered system. If an 
arbitrary ancilla system in a Hilbert space iJ" is introduced, 
the composite state of the principal system and ancilla 
pertains to the Hilbert space H ® H"^. In this case, assuming 
a trivial dynamics on the ancillary, the time evolution of the 
total system is given by p''" (<) = (A(t) (g)/)/9''°(0), where / is 
the identity operator acting on the state space of the ancillary. 
The amount of total correlations in a bipartite state p''" can be 
quantified through the quantum mutual information 

I{pn:=S{p^) + S{p^)-S{D, (1) 

where = tr^p"" and p" = trap*° represent the reduced 
density operators of the system and the ancilla, respectively. 
S{p) = —trplogj p is the von Neumann entropy. Exploit- 
ing the fact that quantum mutual information is a monoton- 
ically decreasing function of time for a Markovian process, 
LFS have proposed a new quantity for measuring the non- 
Markovianity of the dynamical process A(t) from an infor- 
mation theoretical perspective: 

AA(A) = sup / ±np^-(t))dt, (2) 

p=''(0) J(rf/dt)/(p=°(t))>0 0.1 

where the sup is taken over all possible initial states p*''(0). 
Even though this measure has a significant meaning for the 
quantification of non-Markovianity, its evaluation is hard to 
perform due to the potentially complex optimization problem. 

As described in Eq. (2), a dynamical quantum process is 
said to be non-Markovian if ^/(p*''(t)) > at any instant 
of time. Note that the ancilla, unlike the system, does not 
interact with the environment. In other words, the state of the 
ancilla is time independent. Therefore, the time derivative of 
the quantum mutual information can be written as 

= J^s{p'{t)) - pint))- (3) 

For a zero temperature environment, an interesting result can 
be obtained from this equation. Since we take p'^'^(O) as a pure 
state, and the environment starts in the state p'^(O) = |0) (0| for 
the zero temperature case, the total quantum state composed 
of the system, the ancilla and the environment is a pure state 
at any time, leading to 5'(p'"' (i)) = S{p'^{t)). Consequently, 
following the LFS measure, we obtain a non-Markovianity 



criterion without the need of an ancilla: a quantum process 
is non-Markovian if and only if 

Pip^) > pimy (4) 

This condition links the non-Markovianity measure to the 
rate of change of the system and the environment entropies. 
Nonetheless, some peculiar aspects should be noted: the en- 
vironment is initially in a pure state by assumption, but no 
restrictions were imposed to the system. If the system is ini- 
tially in a pure state as well, we have p'*^ pure and therefore 
S{p^{t)) = S{p'^{t)). Such a result does not mean that the 
process is actually Markovian because we need to maximize 
over all possible initial conditions of the system to be able 
to determine the degree of non-Markovianity. An equivalent 
equation for Eq. (2) can be deduced without the necessity of 
an ancilla: 

Mt=ok{^) = sup / —ASse{t)dt, (5) 

^"(0) J{d/dt)ASsa{t)>0 Ot 

with ASse{t) = Sif^it)) - Sip^it)). The advantage of the 
above equation over Eq. (2) is that the maximization is just 
over the initial conditions of the system instead of the initial 
conditions of the composite state of the system and the ancilla, 
that is required to calculate Af{A). It is important to empha- 
size that, for a non-Markovian process, even though the en- 
tropies of system and environment can increase or decrease, 
the main point is connected to the time derivatives of the en- 
tropies. Furthermore, to calculate S{p'^{t)), we do not need to 
worry about the state of the environment. The idea is to max- 
imize AfT=oK over all possible initial system states and, for 
each choice, purify it including an extra subsystem. Because 
the environment is set in a pure state at < = 0, the entropy of 
the system plus the purifier subsystem is equal to the entropy 
of the environment at any time. 

On the other hand, LFS presented a significant simplifica- 
tion for Eq. (2) in ref. [9]. Assuming that H'^ — H and 
p'*''(0) = |^'}(^| where Y^) is any maximally entangled pure 
state of the system and the ancilla, they obtain an easily com- 
putable measure of non-Markovianity: 

AAo(A)- / p{p'''{t))dt, (6) 

J(d/dt)i(p-'-(t))>Q «t 

with p'*"(i) = (A(i) (g) /)|\I')(\&|. In this work, we explic- 
itly show that A/o(A), despite its utility as a witness for non- 
Markovianity, may be misleading and give an inaccurate con- 
clusion about the degree of non-Markovianity of a quantum 
process. Furthermore, we demonstrate that AA(A) does not 
depend on the amount of entanglement shared between the 
system and ancilla because two distinct initial states with the 
same degree of entanglement can give different results. Actu- 
ally, the optimal state is not maximally entangled in general. 
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III. OPEN SYSTEM DYNAMICS 



r{t) is 



A. Phase Damping Channel 

We consider a spin-boson type Hamiltonian Hipo that de- 
scribes the interaction between a qubit and a bosonic envkon- 
ment: 



HiPD ~ 



k k 



<T:,{gkal+ glttk), (7) 



where the first and the second terms are responsible for the 
free evolution of the qubit and the environment, respectively. 
The third term accounts for the interaction between the qubit 
and its environment. We first note that [H, <Tz] =0, which im- 
mediately implies the absence of transitions between different 
energy levels. Thus, the population terms in the density matrix 
of the system are conserved quantities. Here, wq is the tran- 
sition frequency of the qubit and Uk is the field frequency of 
the fc-th environmental field mode. The constant controls 
the strength of the coupling between the qubit and each field 
mode of the environment. While the qubit operator is given 
by the usual Pauli az matrix, the creation operator Ofc and the 
annihilation operator a[., satisfying the bosonic commutation 

relations [afc,a[.,] ~ Sk^k', represent the environment. It is 
worth to stress that this qubit plus environment model admits 
an exact solution. We assume that the composite state of the 
qubit and the environment are initially factorized, that is, there 
exists no correlations between the system and the environment 
at t = 0; furthermore, the environment is initially in its vac- 
uum state p'^{0) = |0)(0| at zero temperature. We consider 
a sufficiently large environment for which we can replace the 
sum over the discrete coupling constants with an integral over 
a continuous distribution of frequencies of the environmental 
modes, i.e., |gfcp dujj{uj). In addition, we sup- 

pose that the spectral density of the environmental modes is 
Ohmic-like, that is 



(8) 



with uJc being the cut-off frequency and 77 a dimensionless 
coupling constant. Depending on the parameter s, the spectral 
density is called subohmic (s < 1), ohmic (s = 1) or super- 
ohmic (s > 1). Under these conditions, the dynamics of a 
single qubit can be obtained in the operator-sum representa- 
tion as 



1=1 

where the Kraus operators Kj (t) are given by 



(9) 



(10) 



with X]i=i = I for all values of t, where / de- 

notes the 2x2 identity matrix. Here, the dephasing parameter 



"(t) = exp 



-/{t')dt' 



(11) 



where the dephasing rate takes the form 

7(t) = rjUJcil + iuJct)'^)~'^^T{s) sin(s arctan(wct)), (12) 
with r(s) being the Euler gamma function. 

B. Amplitude Damping Channel 

In order to discuss the relaxation process, we consider the 
following model Hamiltonian Hjad 

Hi AD = ojQ(J+(T- + ^ ojkalttk + {(t+B + cr-B^), (13) 
fe 

where B = J^k 9kO-k with gk being the coupHng constant. 
The first two terms of Eq.(13) describe the free evolution of 
the qubit and the environment, respectively, while the third 
term accounts for the interaction between the qubit and the 
environment. The transition frequency of the qubit is ujq, and 
(T± denotes the raising and lowering operators of the qubit. 
The index k is used to label the different environmental field 
modes with frequencies cj^, which are mathematically de- 
scribed by the annihilation and creation operators given by Uk 
and aj., respectively. Restricting ourselves to the case of a sin- 
gle excitation, the modes of the environment can be described 
by an effective spectral density of the form 
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27r {luo - uj)^ + A2 ' 



(14) 



where A defines the spectral width of the coupling and is also 
connected to the correlation time of the environment tb by 
the relation tb ~ 1/A. 70 is the time scale over which the 
state of the system changes hy Tp w l/7o- For this form of 
a spectral density, it is not hard to distinguish the weak and 
the strong coupling regimes. The case tr > 2tb corresponds 
to the weak coupling regime where the decoherence process 
is Markovian since the relaxation time is greater than the cor- 
relation time of the environment. On the other hand, the case 
Tp < 2tb corresponds to the strong coupling regime where 
the non-Markovian nature of the environment becomes evi- 
dent. We note that at zero temperature this Hamiltonian with 
the considered spectral density (known as the damped Jaynes- 
Cummings model in the literature) represents one of the few 
exactly solvable models for open quantum systems. In the 
strong coupling regime, the time evolution of a single qubit 
can be expressed in the operator-sum representation as 



p(t)=^M,(i)p(0)M/(i), 



(15) 



where the corresponding Kraus operators Mj (t) are given by 



Mi{t) 











, (16) 



4 



satisfying the condition X]?=i Ml{t)Mi{t) = I for all values 
of t. The damping parameter reads 



p{t) = e 



-At 



n 2 



(17) 



with d = ^27oA - A^. 



IV. NON-MARKOVIANITY FOR A SINGLE QUBIT 

Before starting to elucidate the scaling properties of the 
LFS measure for multipartite systems at zero temperature, we 
consider the case of a single qubit system. In this case, the 
maximization in Eq. (5) can be numerically evaluated because 
the general form of the density matrix ^"(0) depends just on 
three real variables. Explicitly, 



Pii(O) 5ft[pi2(0)] + Q[pi2(0)] 

5R[pi2(0)] -3[pi2(0)] l-pii(O) 



For the dephasing process, our numerical analysis show that 
the maximum present in Eq. (2) is reached for a maximally 
mixed initial state. In other words, the optimal state of the 
composite system (0) is maximally entangled, justifying 
the simpUfication proposed by LFS, as described in Eq. (6). 
However, contrarily to what is expected, the maximum value 
in Eq. (2) for the relaxation process is obtained for a diag- 
onal initial state whose system plus ancilla is not maximally 
entangled. 

As suggested by our numerical investigation, we first take 
the off-diagonal elements of the density matrix to be zero, that 
is, pi2(0) = P2i(0) = 0. In Fig. (1 -a) we plot the possible 
values of the degree of non-markovianity quantified by the 
LFS measure A/'(A), for the PD channel with s = 3, Wc = 1 
and 7/ = 2 as a function of the density matrix population 
pii(O). As can be seen from the figure (1-a), the LFS mea- 
sure corresponds to having jOii(O) — 0.5, implying that the 
optimal initial composite state is maximally entangled. On the 
other hand. Fig. (1-b) presents the results of the same analysis 
performed for the AD channel with 70 = 1 and A = 0.1. It is 
shown that the LFS measure is obtained for pii(O) « 0.4, 
which means that the optimal composite state is not maxi- 
mally entangled in this case. In particular, such a result points 
out that the LFS measure does not generally depend on the 
initial entanglement between the system and the ancilla. 

Another interesting point is related with the dependence of 
the optimal initial state of the system on the parameters of the 
considered environmental model. In Fig. (2-a) and Fig. (2- 
b), we plot the density matrix element pii(O) of the optimal 
initial state for the dephasing and the relaxation processes, re- 
spectively, as a function of the bath parameters Wc and A. The 
results of this analysis demonstrate an important difference 
between the two cases. Whereas the optimal initial state for 
the PD channel does not depend on the bath parameter ujc, the 
AD channel is highly sensitive to the bath parameter A. In 
fact, we have found that as the parameter A gets smaller, the 
optimal state tends to be a maximally mixed one. 
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no. 1 : Non-Markovianity for one qubit as a function of the density 
matrix element pii(O) for the dephasing process (a) with s = 3, 
Wc = 1 and r; = 2, and for the relaxation process (b) with 70 = 1 
and A = 0.1. 
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FIG. 2: The density matrix element pii(O) of the optimal initial state 
of a single qubit system as a function of the bath parameters for the 
dephasing process (a) with s = 3 and r] = 2, and for the relaxation 
process (b). 



V. SCALING OF NON-MARKOVIANITY 

A. Independent Environments 

Having collected all the required tools, we are in a position 
to evaluate the LFS measure for the multipartite case for in- 
dependent environments. First of all, we consider a quantum 
state which is composed of a system of n qubits and an ancilla 
that purifies the system state. Assuming that only n qubits are 
subjected to the environmental noise, and the ancillary system 
evolves freely, the dynamics of the composite system p'^"' (t) 
can be obtained as 

p""{t)^ (£^f"®^)p'"(o)(^^f"®/)^ (18) 

i={l,2} 
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where Ei are the Kraus operators describing AD or PD chan- 
nels for a single qubit, n is the number of qubits, / denotes the 
identity matrix with dimensions of the ancillary system, and 
the summation over the index i = {1,2} runs over all possible 
permutations of the Kraus operators Ef". 
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FIG. 3: Non-Markovianity A/'(A) for the dephasing process (a) with 
s = 3, Wc = 1 and r) = 2, and for the relaxation process (b) with 
7o = 1 and A = 0.1, as a function of the number of qubits. The 
graphs are for independent environments. 

Performing the optimization required for the calculation of 
the LFS measure A/'( A) becomes a very difficult task for a sys- 
tem of two or more qubits due to the significantly increasing 
number of variables involved in these cases. To overcome this 
difficulty, we limit our study to diagonal product initial states 
of the form (O'*(0)®" for multipartite systems. Indeed, this is a 
reasonable choice given the results acquired for the case of a 
single qubit and, moreover, using such diagonal states enables 
us to obtain a lower bound to the degree of non-Markovianity 
of a quantum process. While we plot the LFS measure as a 
function of the number of qubits for superohmic dephasing 
environments with s = 3, = 1 and ?/ = 2 in Fig. (3-a), 
similar results are displayed for the relaxation process with 
7o = 1 and A = 0.1 in Fig. (3-b). These plots indicate a 
linear increase in the degree of non-Markovianity for the con- 
sidered initial states, which proves that, for independent PD 
and AD channels, the LFS measure is at least additive, i.e., 

AfiA®"") > n[M{A)]. (19) 

This peculiar behavior of the LFS measure, despite non- 
intuitive, is simple to be understood in this context because 
both the system and the environment entropies are additive 
quantities when independent environments are considered. 



B. Common Environment 

In this section, we turn our attention to the behavior of non- 
Markovianity when the system is globally interacting with 
common dephasing and relaxation environments. Due to the 



difficulty of calculating the dynamics for many qubits consid- 
ering their interaction with common environments, we restrict 
our analysis to only two qubits. On the other hand, the anal- 
ysis of the two-qubit case is still interesting to infer the main 
characteristics of the degree of non-Markovianity as a func- 
tion of the system scale. We first consider the case where both 
qubits are coupled to a common dephasing bath, described by 
the following Hamiltonian 

HCPD = i5]4"Vi")+^c.,.4afe+^ai")(g,.4+,g;,afe), 

n k n,k 

(20) 

where the index n denotes the terms related to the first {n = 1) 
and second (n = 2) qubit. As a simplification, we only con- 
sider diagonal two-qubit initial states and evaluate the degree 
of non-Markovianity measured by the LFS measure. The re- 
sults of our analysis suggest that the LFS measure N{h) suf- 
fers a very significant decay as compared to the single qubit 
case, having a value of the order of 10~^. 

Next, we focus on the scenario where a system of two qubits 
globally interact with a common relaxation environment: 

HcAD - ^+E^^4«fc+E('^+ '^+^-"'^^)- 

n k n 

(21) 

In this case our investigation reveals that, assuming diagonal 
initial states, the degree of non-Markovianity is significantly 
amplified when compared to the single qubit case, approxi- 
mately turning out to be A/'(A) k, 6.21. Interestingly, among 
the two-qubit diagonal states that we have considered, the op- 
timal one is always the maximally mixed state, independently 
of the bath parameters. Such a finding is rather surprising 
because the optimal initial state is strongly dependent on the 
parameters of the environment for the single qubit case. 

VI. COMPARISON OF A/'( A) AND M) ( A) 

We now compare the results obtained for the LFS measure 
7V(A) and its simplified version A/o(A) for multipartite sys- 
tems, considering independent environments. While the for- 
mer involves a difficult maximization over all possible initial 
states, the latter can be directly calculated choosing a specific 
initial state, which we choose as a GHZ type state. Due to this 
restriction, it is clear that 7\/o(A) underestimates the degree 
of non-Markovianity, and consequently A/o( A) < 7\A(A). In 
Fig. (4-a) we plot the logarithm of the simplified LFS mea- 
sure, ln(A/'o(A)), as a function of the number of qubits for the 
dephasing process with s = 3, Wc = 1 and jy = 2. It can 
be observed that the degree of non-Markovianity measured by 
7Vo(A) decays exponentially as a function of the number of 
qubits. Furthermore, we see that even for very small systems 
(two qubits), the degree of non-Markovianity diminishes very 
significantly. In Fig. (4-b) we make the same analysis for 
the relaxation process considering the parameters 70 = 1 and 
A = 0.1. Our findings demonstrate that, unlike in the case of 
dephasing, the simplified LFS measure A/o(A) might increase 
for the relaxation process. As a result, comparing Fig. (4-a) 
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and Fig. (4-b) to Fig. (3-a) and Fig. (3-b), we conclude that 
Afo{A) might be a misleading quantity for determining the de- 
gree of non-Markovianity, despite the fact that it is an easily 
computable witness of non-Markovianity. 
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FIG. 4: Logarithm of the non-Markovianity ln(A/'o(A)) for the de- 
phasing process (a) with s = 3, Wc — 1 and rj — 2, and the non- 
Markovianity A/'o(A) for the relaxation process (b) with 70 = 1 and 
A = 0.1, as a function of number of qubits. The graphs are for inde- 
pendent environments. 



cient method for calculating the degree of non-Markovianity 
due to the fact that the Hilbert space, where the maximization 
is evaluated, does not include an additional ancillary system. 
We provide an extensive analysis of the non-Markovianity 
for a single qubit and determine the optimal initial states 
of the system required for the evaluation of the LFS mea- 
sure, as a function of the parameters of the environment. We 
demonstrate that the degree of non-Markovianity might in- 
deed increase with the number of the qubits in the system 
for independent environmental interactions and moreover ob- 
tain a lower bound to the LFS measure in this case, namely 
A/'(A'^") > n[M{A)]. Furthermore, we investigate the non- 
Markovian behavior of two qubits interacting with a common 
environment and show that the amount of non-Markovianity 
is amplified as compared to a single qubit for the relaxation 
process. Finally, we present a comparison between the LFS 
measure Af{A) and its simplified version Afo{A) for indepen- 
dent environments. Our examination clearly suggests that al- 
though A/'o(A) can be used as a non-Markovianity witness, it 
may give incorrect results for multipartite systems as a non- 
Markovianity measure due to the lack of optimization in its 
evaluation. 
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VII. SUMMARY 

We study the non-Markovianity of independent and com- 
mon dephasing and relaxation processes for multiparitite 
quantum systems. Considering zero temperature environ- 
ments, we show that no ancillary system is required to quan- 
tify the degree of non-Markovianity of a dynamical quantum 
process since the LFS measure J^{A) can be calculated by 
the difference of the time derivatives of the system and the 
environment entropies. This simplification results in an effi- 
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